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We calculate the single logarithmic contributions to the quark singlet and gluon matrix of timelike 
splitting functions at all orders in the modified minimal-subtraction (MS) scheme. We fix two of 
the degrees of freedom of this matrix from the analogous results in the massive-gluon regularization 
scheme by using the relation between that scheme and the MS scheme. We determine this scheme 
transformation from the double logarithmic contributions to the timelike splitting functions and 
the coefficient functions of inclusive particle production in e + e _ annihilation now available in both 
schemes. This approach fixes two of the four degrees of freedom, and a third degree of freedom is 
fixed by reasonable physical assumptions. The results agree with the fixed-order results at next-to- 
next-to-leading order in the literature. 

PACS numbers: 12.38.Cy,12.39.St,13.66.Bc,13.87.Fh 

I. INTRODUCTION 

According to the factorization theorem, the cross section T (x, Q 2 ) = Q 2 (da h /dx)(x,Q 2 ) for the inclusive produc- 
tion of a hadron h carrying a fraction x of the available energy in a process with an energy scale Q much greater than 
the asymptotic scale parameter Aqcd of QCD can be calculated by the convolution 



D h J-,^ f ), (1) 



z 



where F a (z, fJ- 2 /Q 2 , a s ([i 2 )j = Q 2 (da a /dz) (z,Q 2 ,fJ/j\ is the equivalent process-dependent factorized cross section 

for the production of a parton a carrying away a fraction z of the available energy, which contains all subprocesses 

with energy scale greater than the arbitrary factorization scale fXf, (z, fJb 2 A is the factorized fragmentation function 

(FF) for the fragmentation of a parton a to a hadron h carrying away a fraction z of the energy of this parton, which 
contains all subprocesses with energy scale less than fif, and a s = a s / (2ir), with a s being the strong-coupling constant. 
These partonic cross sections are perturbatively calculable. I.e. the series 

Fa (z,^,a.(fi)j = £ aJfciJJFW (z,^j (2) 

approximates F a (z,^ 2 /Q 2 ,a s (jJ, 2 )\ for sufficiently small values of a s . We refer to this approach to calculations, 

namely expanding in a s with coefficients that are independent of a s , as the fixed-order (FO) approach. The scale 
Q will be loosely defined to be the scale which fj,f should be chosen to have the same order of magnitude as in 

order that the pcrturbative scries for F a (z, (Jj/Q 2 , a s ([i 2 )j be as convergent as possible. To be explicit, this is 

because the coefficients F^ (^z,n 2 /Q 2 ^J in Eq. ([2]) grow like ln n+p (p 2 /Q 2 ) as [J. 2 /Q 2 — > oo, where p is an integer 

that is independent of n. The predictive power of this approach follows from the facts that the FFs are process 
independent and the dependence of the FFs on fj,f obeys the Dokshitzer-Gribov-Lipatov-Altarelly-Parisi (DGLAP) 
evolution equation [l|-[4| , 



din- 2 



^ ! M 2 / ) = E/ 1 f^(^^/))^(^ 2 /) 



(3) 
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where P a p(z,a s ) are the a — > splitting functions, which are perturbatively calculable in the FO approach for 
sufficiently small values of a s , the perturbative series taking the form 

oo 

P a p(z,a s ) = J2<PiT 1) ( z )- ( 4 ) 

n=l 

However, the FO approach fails when x is too small, due to the presence of unresummcd large soft-gluon logarithms 
(SGLs) in the timelike splitting functions and in the hard partonic cross sections. This means that small- a; measure- 
ments cannot be used to provide constraints on FFs at small values of z. They also cannot be used to improve the 

FFs at higher values of z because, according to Eq. ([T]), the cross section at x depends on the FFs D 1 ^ (^z,/i 2 ^j at all 

z values in the range x < z < 1. 

To improve the accuracy at small values of x, the SGLs of each class appearing in the FO expressions must be 
determined to all orders. The double logarithms (DLs), being the largest SGLs, arc known to all orders in the 
MS scheme for the timelike splitting functions Q and the coefficient functions for inclusive hadron production in 
e + e~ annihilation @, 0]- The single logarithms (SLs) in the splitting functions to all orders are known only in the 
massive-gluon regularization (MG)scheme @ . Because the FO approach and the resummed SGLs can be consistently 
combined as discussed in Refs. [H,!! E3] to give an approach which can describe data from the smallest to the largest 
values of x and because FO calculations in the MS scheme are known to next-to-leading order (NLO) and beyond, 
while those in the MG scheme are usually not, it is necessary to determine the SLs in the MS splitting functions. 
Furthermore, FFs are usually determined in the MS scheme. 

In this paper, we first partially constrain the complete SL contributions to the MS splitting functions using three key 
ingredients: firstly, the DL contributions to the splitting functions in these two schemes; secondly, the SL contributions 
to the MG splitting functions; and thirdly, the DL contribution to the scheme change between the MG and MS schemes. 
The third ingredient can be obtained because the DL contribution to the gluon coefficient function of e + e~ annihilation 
is known in the MG scheme and we recently calculated the same quantity in the MS scheme @, 0] ■ To completely 
constrain the SL contributions to the MS splitting functions, we then introduce some reasonable assumptions that fix 
the next-to-lowest order of the scheme change: We demand that our results are consistent with the next-to-next-to- 
leading-order (NNLO) splitting functions, and also that the matrix exhibits certain non-singular properties at small 
values of z. 

This paper is organised as follows. In section ILT1 we discuss the calculations of factorized cross sections in general. 
In section 11111 we introduce SGLs and present the DLs in the MG and MS schemes for the coefficient functions in 
e + e~ annihilation and for the timelike splitting functions. We formalize the relation between calculations in different 
schemes in section UVl In section [V] we use these results together with the SLs in the timelike splitting functions in 
the MG scheme determined in Ref. Q to determine the SLs for combinations of the splitting functions in the MS 
scheme. Finally, we present our conclusions in section IVIII 



II. GENERAL FACTORIZED INCLUSIVE HADRON PRODUCTION CROSS SECTIONS 



In this section, we consider the general structure of the calculations of factorized cross sections that will be useful 
later. We will find it convenient to work in Mcllin space, defined by the (invcrtiblc) Mcllin transform 

/(w)= / dxx»f(x), (5) 
Jo 

because x-space convolutions reduce to simple products. In particular, removing the superscript h from now on, Eqs. 
(P), ([3]), and (0| respectively become 

T{u,Q 2 ) = Y t Fa(^,^,a.(jt)^D a {w, l ? f ), (6) 

-^-^D a (u,fi 2 f ) = Y^P*l}(u,a s ((J? f ))Dp((j,fj, 2 f ) , (7) 
nfi f p 

OO 

P a p(u;,a s ) = E a " P S _1) H- ( g ) 

n=l 

According to Eq. ([6]), the cross section is invariant under any change of parton basis F a F a = J2p F/3(Y~ 1 )p a 
and D a — > D a = Y a pDp, where Y is any invertible matrix which is independent of ui, p?p and Q 2 . In matrix 



3 



notation, F = FY^ 1 and D = YD. For example, the SU(n/) symmetry of the DGLAP equation in the MS scheme for 
Uf active flavours of quarks and the charge conjugation symmetry of QCD imply that P is reduced to block-diagonal 
form when the parton basis is chosen such that the FFs consist of the quark singlet component, 

, "/ 

Dv = —J2(D qj +D qj ), (9) 

Uf .7=1 

with qj (q,j) being the (anti)quark of flavour J, the quark non-singlet component, 

D qj ,vs = D qj + D qj - D s , (10) 
the valence-quark singlet and non-singlet components, and the gluon component, D g . In this basis, for 



in Eq. , we have the 2x2 matrix 



while, for D = D qj ^s, we have the single flavour- independent quantity P = Pns> an d simlarly for the valence-quark 
singlet and non-singlets. 

An alternative basis, which is used in some applications and will be needed later, is that in which the LO splitting 
function matrix is diagonal, i.e. 



and 



where, defining the projectors a, /3, and e by [111 ] 



p(0) p(0) p(0) p(0) 



p(0) p (0) ' ^ p (0) p (0) ' p (0) p (0) 



D+ = (l- a)D^ - PD 9 , 



(15) 



(16) 



we have 



and, for all k > 0, 

P W = apW+/3P«+cpW + (l-a)J^), 
p (fe) = p (fe) _ Ap(fe) 1 - a p (k) 

H I SS Eg 

p(fe) _ p(Jfe) , p (fc) _ p (fc) 

pM = pW _ (p<M _ 2 P W) = pW ( P W _ 2 P «) . (17) 

Note, of course, that P±l = by definition. In one important simplification of QCD, namely Af = 4 super Yang-Mills 
theory, this basis is actually more natural than the basis of quark singlet and gluon because the diagonal splitting 

(k) 

functions P±± can be expressed in all orders of perturbation theory as one universal function with shifted arguments 
0. 

In general, because both the Mellin transform and the change of parton basis are invertible, we will not specify 
whether the x-space convolution of two x-space functions or the product of their Mellin transforms is being calculated, 
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nor which parton basis is being used, nor whether only a subspace of the full parton space (achieved by setting 
combinations of FFs to zero) is being considered, but simply write Eqs. ((T|) and © as 

T = FD, (18) 

and Eqs. © and © as 

d 



-D = PD. (19) 



Inclusive particle production in e + e annihilation provides a simple example of this formalism. In this case, Q is 
conveniently chosen to be the cm. energy, and the cross section takes the form 

n f 

F =Y. F u,^D qj ,ns + F s Dj: + F g D g , (20) 



J=i 



where 



F ^,ns [u,^,a a (ji 2 f ^j = Q 2 a (Q 2 )N c Q qj (Q 2 )C m ^o,^,a s Qi 2 ) 
F E L,^,a s (/4)j = Q 2 a (Q 2 )N c n f (Q(Q 2 ))Cx L,^,a s (fi) 



F g L,^,a a ^ 2 f )\ = Q 2 a (Q 2 )N c n f (Q(Q 2 ))C g L^,a s (fx})\, (21) 

with <Jo(Q 2 ) being the lowest-order (LO) cross section for the process e + e~ — > N c the number of quark colours 

in QCD, Q qj (Q 2 ) the effective electroweak charge of quark qj, and (Q(Q 2 )) = Y^jLi Qq.i(Q 2 )/ n f- Note that Eq. (f2"Tj) 
is, strictly speaking, dependent on Q through Al^/Q 2 (in Q qj {Q 2 )) as well as through [i 2 /Q 2 , but these dependences 
are not shown for brevity. The coefficient functions Cx (X = NS,S,g) in the FO approach in Mellin space may be 
found, e.g., in Ref. [13| . It will be convenient later to write 

F = Q 2 a {Q 2 )N c n f (Q(Q 2 ))CD. (22) 

For example, for the quark singlet and gluon contribution in Eq. (|20[) . 

F = F S ^ E + F g D g , (23) 

D is given by Eq. (|11[) . and 

C = (C^C g ). (24) 
We will set [if = Q for simplicity, in which case it is convenient to define 

C x {uj,a s ) = C x {uJ,l,a s ) (X = NS, £,<?). (25) 

III. SOFT-GLUON LOGARITHMS 

Since the non-singlet inclusive partonic production cross sections i^ gj) NS and the non-singlet splitting functions 
are free of SGLs, they do not concern us, and so we will not discuss them further. From now on, inclusive particle 
production cross sections will be assumed to take the form in Eq. (|23p . The inclusive partonic production cross 
sections F calculated in the FO approach may exhibit a singular behaviour in Mellin space as uj — > 0. This is caused 
by SGLs, which grow like \/uj p for p > 1. In a; space, these SGLs take the form of quantities that grow like ln p_1 x 
as x — > 0. Such strong singularities are non-physical and become weaker or even disappear after being resummed to 
all orders. The resummed SGLs in F take the form of the series 



oo 

f=j2 (—\ 

m=0 
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For such a series to converge, at least asymptotically, it is necessary that a s «l and u> — O(yfal). The DLs, namely 
those SGLs for which m = in Eq. (|26p . of the inclusive partonic production cross sections for e + e~ annihilation in 
the MS scheme, when D is given by Eq. (fTTj). take the form 



C = (l,C g DL ), 



where 



with 



2C F 



4 7 (cj, a s ) + us 



7(w, a s ) = + V^ 2 + 16CUa s ), 



(27) 



(28) 



(29) 



and Cns = 1- They were also determined in Ref. |14j in the MG scheme, indicated in this paper by an overline, to be 

(30) 



C=(1,C° L ), 



where 



C g (w,a s ) = — 



47 (w, o«) +a; 



(31) 



and Cns = 1- 

The rcsummation of the SGLs in P take the form of the series 



'-EG)"*" 1 (5 



(32) 



The full DL contribution to P, namely the SGLs for which m = 1 in Eq. (|32|) . will be written as P DL = 
(a s /o;)p[ ](a s /w 2 ). When L> is given by Eq. it is given in the MS scheme by 



where 7 is given in Eq. (|29p and 



P UL (w,a s ) = A 7 (w,a s ), 





1 



(33) 



(34) 



-DL 



which obeys the projection operator property A 2 = A. For the quark non-singlets, P =0. The DLs in P in the 
MG scheme are the same as those in the MS scheme, i.e., when D is given by Eq. (fTTj) . 



-DL 



P (u,a s ) = Ay(uj,a s ), 



(35) 



-DL 



and P = for the quark non-singlets. 



IV. GENERAL SCHEME CHANGES 



Results in one scheme, such as the splitting functions in the MS scheme, may be obtained from the analogous results 
in another scheme, such as the MG scheme, once the relation between the two schemes is known to the appropriate 
accuracy. To obtain the form of this relation, let F (D) and F (D) be respectively the partonic cross sections (FFs) in 
any two different schemes. Since the cross section is scheme independent, then as well as Eq. (fT8|) we have T = FD. 
Comparing this last result with Eq. (fTSf and treating D a (z,iij) as arbitrary functions, we find that 



F 
D 



FZ,_ 



(36) 
(37) 
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where, in Mellin space, Z is an invcrtible matrix that depends on u> and fij. Note, therefore, that Eqs. ([36]) and ([37)) 
are generalizations of the change of parton basis considered just after Eq. ([7]). The (matrix of) splitting function(s) 
P is defined to be that which appears in the DGLAP equation in the new scheme, which emerges from Eq. (fT9]) by 
substituting D and P with D and P, respectively. Then, it follows from Eqs. (|T5|) and (|3"T|) that the relation between 
the splitting functions in two different schemes is given by 

P = Z-'PZ - Z- 1 -^—. (38) 
a In 

Now, consider a general expansion of perturbativcly calculable quantities, such as F a and P, in some variable 
x{w, a s ), with coefficients that depend on y(uj,a s ), i.e. 



F = J^x n F^(y), 

n=0 
oo 

p = Y, xnp{n ~ 1} (y)- ( 39 ) 

n = l 

For example, in the FO approach (Eqs. ([2]) and ([8])), x — a s and y = oj in Eq. (|39|) . while in the SGL approach (Eqs. 
(|26l) and (J32J) ) , x = a s /uj and y = a s /uj 2 in Eq. (|39|) (excluding terms that are non-singular as ui — > 0). We restrict 
our schemes to be such that, if the perturbative series for F begins at 0(x n ), the perturbative series for F also begins 
at 0(x n ). Thus, 

OO 

Z(u,a s ) = J2 z{n} (y> n - (40) 

n=0 

Note that, in Eq. ([38]), the first term Z^PZ = 0(x) while the second term Z~ 1 dZ/dln fj? f = 0(x 2 ). Thus 

pM = ^{o}-ip{°} z {o} Thc rcsult p(0) _ p(°) nQ i ongcr hc-ids in general, but rather if and only if Z^ commutes 

with P^°\ This is trivially the case in the FO approach because the schemes used in the literature are (usually) such 
that Zl°y(y) = 1, i.e. 

OO 

Z(u,a s ) = l + J2 z(n) ("K- (41) 
n=l 

However, in the SGL approach, where 

m=0 

we must allow for the possibility that Z^(y) is any function of y. We will see later that Z^ does in fact commute 
with P , at least for the MG and MS schemes. 



SINGLE LOGARITHMS IN THE SPLITTING FUNCTIONS 



Thc SL contributions to the timelike splitting functions have already been calculated in thc MG scheme [8j, [1J, [15| 
and are given by 



pSL 
* EE 

-^E S 



^SE 



5 SL 
99 



o, 

2C F 
C a 



—ST 1 

P 99 + Q(nC A + 4n f T R )a s 



1 Ih/Tr 2C F n f T R 



6 3 C A 



3 



-T R rifa s , 

lLu 3 (nC A + ^n f T R ) 2u[55C A (2-/ + uj)+ACAnfT R (6j + 5uj) + 8C F nfT R io] 2 
~6 (4 7 + w) 3 as ~ 3 (2 7 + w)(4 7 + cj) 3 ° s 

16 (HCj + 12Cln f T R + leCACj.n/Tfl) 3 
3 (2 7 + w)(4 7 + W )3 0s ' 



(43) 

3C F a s , (44) 
(45) 



(46) 
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where 7 = j(oj,a s ) is given by Eq. (|29l . In Eq. (|44|) . we have taken the opportunity to correct some obvious 
typographical errors in Eq. (38) of Rcf. @ 35|. It is the goal of this section to perform the scheme change given in 
Eq. (|38|) on Eqs. (|43| -(j46 |) in order to constrain and then to attempt to determine the SL contributions in the MS 
scheme. 

We first calculate Z^ from the DLs in the coefficient and splitting functions: with the help of Eqs. (|27|i and (|30l) . 
Eq. (f3"r?)) becomes 



(l,C g DL ) = (l,C7) 



7 [o] 7 [o] 

7 \0] 

^gY, ^99 




(47) 



Using this result to eliminate Z^ and gives 



z [0] 1 ( g 'V: ( •/ ~ \ ;.• X 



,[0] 

J 9 s 



[0] 



(48) 



Next, we note that, because dZ/dln/sj is free of DLs as discussed immediately after Eq. (|4H)) [36j . the DLs in Eq. 
obey 



pDL = z [o]-ip^ z [o] 



(49) 



Using Eqs. (|33)l and (|35|) . we find that 



[Zl°\A] = 0, 



(50) 



i.e. 



f(i-OS- 




, z [o] 



2C F 7 [0] 
7 [0] 



2C F y[0] 



[o] 



(51) 



from which wc find that 



z [ f = 
= 

99 



0, 

2C F 



DL 



(52) 



Thus, finally, 



(l 



z n = 



(53) 



We note that 



,[0]-i _ 



1+' 
1+ 



-C DL 



1+s 



(54) 



by taking C^ L O C^ L , as expected. 



is obtained from Z^ 1 

We are now in a position to constrain the SLs in P. Using Eq. 
in the MS splitting functions in Eq. ([38)) are given by 



50)) . the SLs (divided by (a s /ui) 2 for convenience) 



P [i] =zW- 1 P 11] zW+/3 q uj 2 ZM- 



da, 



(55) 
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where (3q = (11/6)C A — (2/3)T R n,f is first coefficient in the series (3{a s ) = — X)^Lo /^" a ™ +2 that determines the running 
of the coupling via da s (fi 2 )/d\nfi 2 = f3(a s (n 2 )) and we have defined 



2CV 1 7 [1] 2C F(7 W yK\ yW (2C F \ 2 ( 1 ,\ 7 [1 

1 _2CV 1 7 [1] 



(56) 



From Eq. (|32p and the definitions that follow it, we have Pjg = (a s /u>) 1 "/(uj,a s ). Explicitly at SL accuracy, Eqs. 
(1551 and (ED read 



-DL 

" g 

"SE - -7T C 4 r DL 
1 + 2C7°g 



pSL _ *-"g *-'g -p SL , nSL 
j w — — — r „r. t -fiss, 

-^DL x 2 



1 1 C4 r»DL / , -DL 7T DL \ (C -C DL , 

, L 9 ~ L 9 J5 Sh TJ SLN 



2C f 9 

pSL _ ^_2CfJ19 "p SL , pSL 
2C f 9 

/~*DL 7^ DL , 1 ! C A /^DL 

pSL _ pSL ° 3 °g -pSL fl 2 d , 1 + 2C F °g R SL 

p gg " p gg + 1+ c^ c dl P ^ + A)a s — In DL + , (57) 

2Cf g 2C F g 

with the definitions P SL = (a s /u;) 2 pM (a s /w 2 ), which are the complete SL contributions to the splitting functions, and 
P SL = (a s /w) 2 pW (a s /u; 2 ). Although our results for the SL contributions to the splitting functions in Eq. ([57|) depend 
on the unknown components of the matrix P SL , it is clear from its form in Eq. (|56p that two of the four degrees of 
freedom of the matrix P SL are completely fixed. For example, these can be taken as any two of P^ — (2Cf /GU)P|e, 
P 9 s g L + (2C F /C A )P^, and the trace P|| + P^. 

We note that, interestingly, R^ = 0(a\) because, with this choice, 

Pih = 32CA zf Rnf aUO(a% 
SL . UC A + 4T R n f 2{llC A + 4C A T R n f -8C F T R n f ) 2 



P 99 



6 s 3w 2 

8(33C^ + l2C A T R n f - 20C A C F T R n f ) 
3u? 



a 3 s + 0(a% (58) 



whose 0(a 2 ) terms are consistent with the results of Ref. and whose O(a^) terms are consistent with the results 

of Refs. [II GJ]. Note also that zf^ = (although we know from Appendix B of Ref. [3 that Z [ ^ ^ 0). We 
will return to this point in the next section, where we investigate the effect of physical constraints on the remaining 
undetermined degrees of freedom on our results. 



VI. PHYSICAL CONSTRAINTS ON THE SPLITTING FUNCTIONS 

In this section, we further constrain the SL contributions to the MS splitting functions by exploiting some pysical 
properties of scheme-dependent quantities in physical schemes such as the MS scheme. 

According to Eq. (|56]l. P|^ is the only component of the splitting function that is so far completely unconstrained, 

even to O(a^), since RJ£ g also depends on the remaining three components of Z^\ and and Zgg are unknown. 
Fortunately, this degree of freedom is fixed by assuming the absence of u — > singularities for all values of the 
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factorization scale in D_ defined in Eq. (|T6|) . By inspection of the DGLAP equation in this basis, it follws that the 
splitting functions P__ and P |_ are found to be free of u> — > singularities, i.e., neglecting all non-singular terms, 



P— = P-+ = 0. 



(59) 



This assumption is expected to be true to all orders. It is certainly true for the DL contributions to the timelike 
splitting functions, for the SL contributions in the MG scheme given in Eqs. ([43|) - (l46l). and through NNLO (l7l. [l8j. 
as we verified in this paper. Moreover, it is true through NNLO in the spacelike case (20l l2i| and holds for the leading 
and next-to-leading singularities to all orders in the framework of Balitski-Fadin-Kuracv-Lipatov (BFKL) dynamics 
j23 - |25| . a fact that has been exploited in various approaches (see, for example, the recent papers (2(| [27j and the 
references cited therein). We note that the timelike splitting functions share a number of simple properties with their 
spacelike counterparts, e.g. the LO splitting functions are the same, and the diagonal splitting functions to all orders 
grow like lnw as u> — > oo. 

Using the relations between the two bases in Eq. (jTTJ) , Eq. ([5TJ1) implies that 



-P 



- € -P 



[19 ■ 



(60) 



where, through the SL level, which is all we need 

2C F ' 



c 
o 



C A 



1 + w - + 



1 1 njT R _ 2 C F n f T R 
6 3 C A 3 C\ 



(61) 



which, at any order k, relates the two most singular terms in the off-diagonal splitting functions P^J and PVg with 

nnH 
EE dnQ r 99 



those in the diagonal splitting functions PyS and Pgg . Using Eq. , the assumption in Eq. ([5^]) implies that 



R 



SL 



flE 2C F 

which is already satisfied by the form in Eq. (|56[) . and 



A Rlk = o, 



(62) 



R 



SL 



^Rf a = \a s C F [C\ + 2n f T R (C A ~ 2C F )} 



=;DL 



C» L (C A C^ + C F ) 



(C A C® L + 2C F )(C A C^ + 2C F ) 



(63) 



which turns out to fix the following combination of components 

"2C 



SL _ 2C F SL [0] 
9 Ca 99 



2C F [i] 
~Ca~ Z ^ 



C A I Z ^ 



C A 99 



(64) 



We can now write Eq. (|56[) in the form 



la s C F [C 2 A + 2n f T R (C A ~2C F ) 











Cf L (C A C° L +C F ) 



(CaC" ij +2C f )(CaC°^+2C f ) 




2C F \^ 
C A J 




2C F 


' Zl°9 ] 


C A 1 



which shows that does not affect the evolution of the combination (2C f /Ca)D 9 — Dt,, reducing the dependence 
of the evolution on the unknown quantity Z g j.. 

Note that Eq. (|59]l (or, equivalently, Eq. (|60|)) implies that the determinant of P vanishes. Since the results in Eq. 
([57t with i?W given by Eq. (|55|) imply that the trace is non-zero, this means that one of the eigenvalues is zero and 
that the other one coincides with the trace. 

In order to complete the check of Eq. (|58|) against the FO results in the literature, we need to consider the remaining 
splitting function, for which we find 



;>SL 

s 9 



—3C F a s + 



12C A C F 2 ,(29C A C F +4C F T R n f (C A -C F )) 3 
~ 16 3^ fls 



0(a 4 s ). 



(66) 
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Here again, the 0{a 2 s ) terms are in agreement with the results of Ref. [16[, while the 0{o? s ) terms are in agreement 
with the results of Ref. [l8| . 

We find that our resummed results exhibit the following u 



behaviour: 



pSL 



pSL 



pSL . 

99 



AC F T R rif a b J 4 4C F rifT R 

3C 3/4 3Ca < 

8C 2 F T R n f a¥ A C F (11C; 



■il P99 zc F 



7/4 



3C A 
2T R n f C A 



'CO 



6C\ 



a* — Z 

3 1 



'] C A 

4C A n f T R - 2AC F n f T R ) {1] pf} ( 2C F \ 



1/4 5/4 

tie 



[o] 



3 y/u 

AC F T R nf a b J 4 



3C" 



3/4 

A 



AC A nfT R - 8C F n f T R 
12C^ ' 



Z 



Pgg 



^99 



(67) 



These limits imply the following nontrivial relation among the SL contributions to the MS splitting functions: 



2C f p sl 
~Ca~ ES 



pSL 



2Cp 
~~Ca~ 



pSL 



2C F p SL 

Ca 99 



(68) 



J, which can be checked using Eqs. (|43|) and 



Equation (|68|) is also obeyed by the SL contributions in the MG scheme [ 
PI?)) , and by the DL contributions, which are the same in both schemes. 

It is interesting to observe that Eq. is also true for all values of to when the choice C A = C F = rif is made, 
which corresponds to an Af = 1 supersymmetric theory. Supersymmetry relations like the one in Eq. (|68[) were first 
introduced in Refs. 0, [28| at one loop, then discussed at two loops in Refs. [3(| and, very recently, at three loops 
in Ref. [l8| . Accidentally, as shown in Ref. @ in the timelike case, Eq. (|68p also reflects the fact that an observable 
like the multiplicity ratio in quark and gluon jets is scheme independent. 

As we have seen, the only undetermined quantity appearing in our formulae, is not constrained by any physical 
conditions. In addition we noted that = 0( a t)- This suggests that is an artefact of the MS scheme. 



VII. CONCLUSIONS 



In this paper, we presented the SL contributions to the MS splitting functions. Two of the degrees of freedom 
in the flavour-singlet matrix were determined from the SL contributions to the MG splitting functions and the DL 
contributions to the e + e~ coefficient functions in both the MS and MG schemes. One of the remaining two degrees 
of freedom was constrained by using certain non-singular properties of the flavour singlet matrix at small values of 
u>, which have been investigated only in the spacelike case so far (see, e.g., Ref. [HI). Nevertheless, both eigenvalues 
are determined analytically in closed form. Our results are in agreement with very recent calculations of the splitting 
functions in the MS scheme at NNLO [13, EH in the FO approach, and also with general physical requirements such 
as supersymmetry. 

Knowledge of the complete SL contributions to the splitting functions formally improves the theoretical description 
of the evolution of FFs at small values of u> and thus facilitates the extraction of FFs from experimental data at small 



values of x in global fits. To date, such global fits have been performed to NLO in the FO approach [31| 



Oi 



calculation of the SL contributions can be incorporated into such fits using the consistent approach of Ref. [5( , which, 
together with the DL contribution to the e + e~ coefficient function determined in Ref. @, allows for a description of 
the experimental data from the largest to the smallest x values. We recall that the NLO splitting functions contain 
also sub-SLs (sSLs), namely the SGLs for which m = 3 in Eq. (|32|) . proportional to aj./u), but the complete sSL 
contributions to the splitting functions are unknown. In the SGL+FO(+F05) scheme defined in Ref. Q, these sSLs 
are, therefore, simply subtracted at this logarithmic order of accuracy. Alternatively, these unresummcd sSLs can be 
replaced by a simple matrix of sSL functions, which are non-singular as u — > 0, but whose FO expansions start with 
the NLO sSLs. An example of such a matrix of functions is a s ~fA/ (2C A ). 

Note added. After the completion of this work, there appeared a preprint [3~ij containing an alternative calculation 
of the SL contribution to the MS splitting functions, with which we found agreement, confirming both our approach 
and the approach of that article. We stress that our approach highlights the relation between the splitting functions 
in the MS scheme and those in the MG scheme, the latter scheme being important in the scheme independent ratio 
of gluon to quark jet rates @. We also showed explicitly that our results are consistent with physical constraints. 
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Finally, we obtained closed forms for the splitting functions (up to Rg%) which we explicitly used to verify the physical 
result in Eq. ((551) . 
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